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Abstract 

We consider the polarization tensor of a color charged gluon field in 
SU(2) in the background of a homogeneous color magnetic field and cal- 
culate its projection onto the lower (tachyonic) state. We obtain a quite 
simple representation of the polarization tensor in terms of double para- 
metric integral. We provide examples for the numerical evaluation of the 
polarization tensor in this case. Also, we consider the corresponding ex- 
pressions at high temperature and calculate the charged gluon magnetic 
mass in this limit. We conclude that radiation corrections stabilize the 
ground state. Applications of the results obtained are discussed. 



1 Introduction 

Color magnetic fields are actively discussed in some areas of modern physics. So, 
the spontaneous generation of such fields prior to the electroweak phase transition 
is considered in [U |2j |3j 0] as candidate source for the intergalactic magnetic fields 
observed today [5j [6] . Strong magnetic fields are considered to be created for a 
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short time in heavy ion collisions at LHC [7J |8]. For instance, their influence on 
the temperature of the de-confinement and chiral phase transitions is discussed. 

Especially interesting quantities are the temperature induced masses which 
determine the characteristics of the quark-gluon plasma. Without magnetic field, 
the Debye mass (defined as m 2 ^ = n 4 4(T, fc 4 = 0, k — > 0) from the gluon 
polarization tensor) is of order mpj ~ gT (where g is the gauge coupling constant 
and T is the temperature) and the magnetic mass, defined in a similar way 
from the transversal components of the polarization tensor, is of order m m ~ 
g 2 T. In this way, both, chromo-electric and chromo-magnetic fields are screened 
in hot QCD. These masses were calculated by different methods, for instance 
analytically from the 1-loop polarization tensor in [9], or by lattice calculations 
in [10]. We remind the corresponding results in QED [11], resulting in ~ eT 
with the electric charge e, and a zero magnetic mass, m m = 0. As a consequence, 
the magnetic fields remains long range in hot QED. 

In the presence of a color magnetic background field, these characteristics 
become modified. In this case it is reasonable to rewrite color gluon field, 
(a = 1,2,3), in terms of charged, = ^(V* ± iV 2 ), and neutral, = A^, 
ones. The Debye mass becomes a function of the magnetic field, ttid = to.d(T = 
0)f(T/y/gB), where f(T/y/gB) is some function turning into unity for B — > 0, 
but reducing mjy for any finite B (see [12] for the neutral and [13J for the charged 
gluons). For the magnetic mass a different picture is observed. For the color 
neutral gluon, in [12] the magnetic mass m m was found to vanish, m m = 0, and 
the account for polarization tensor results merely in a renormalization of the 
external field like in QED. 

It is interesting to note that the structure of the polarization tensor makes the 
characteristics of the gluon field at finite temperature in a magnetic background 
field quite different from those without the field. If one considerers the polariza- 
tion tensor in one-loop approximation at finite temperature without background 
field, a fictitious pole appears [9j HI]. One needs to do re-summation (e.g., the 
hard thermal loop re-summation [15]) to get rid of it. The emerging masses 
are mo ~ gT for the electric and m m ~ g 2 T for the magnetic components. In 
opposite, in a magnetic background field, there is no fictitious pole already at 
the one loop level at any temperature. As shown in [12], the projection of the 
right hand side of the Schwinger- Dyson equation, D^ 1 {k 2 ) = k 2 — II 2 , does not 
change sign. For example, for high T, the projection on the ground state takes 
the form {D^ik 2 )) = k 2 {l + 5.7 T/ y 'gB). Hence, in the magnetic background 
no re-summation is necessary. 

In the present paper we use the representation of the polarization tensor of the 
color charged gluons obtained in [13] and calculate the projection on the ground 
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state. This is the tachyonic state on tree level. This state is given by 



i 


V o j 



o) (i) 



(see eqs. (27) and (32) in [16]). An interesting question to answer here is whether 
the polarization tensor is able to remove the instability of this state. 

For the color charged gluon, the corresponding projection of the polarization 
tensor is still not calculated in a sufficiently explicit form. The present paper is 
devoted to closing this gap. This case requires a special consideration different 
from a general procedure developed already in [13J. Thereby we pay special 
attention to the structure of the ultraviolet divergent pieces and confirm their 
independence of the background field. In addition, the high temperature case 
will be investigated. 

In the next section we adduce the necessary formulas and notations for the 
polarization tensor. The actual calculation for the zero temperature is carried out 
in the appendix. In sect. 3 we consider the high temperature limit for the polar- 
ization tensor and calculate the magnetic mass. The conclusions and discussion 
are given in the last section. 

2 Calculation of the polarization tensor 



We use the representation of the polarization tensor given in [13J. In momentum 
representation, the initial expression reads 

/dk 
Jc^y {^\upG vv '{p - k)Y X 'v>p>G pp >{k) 

+ (p - k) X G( P - k)kyG(k) + k X G( P ~k){p- k)yG(k)} 
+n tadpol ^ (2) 

where the second line results from the ghost contribution and the tadpole contri- 
bution is given by 

nS po1 = / ^ {6xxG pp (k) - G(k) xy } 

+ / 7^4 tovG*,(p) + G ^(P) - 2Gaa'(p)} • (3) 

The vertex factor, 

T Xup = (k - 2p) p 8 Xv + 5 pu {p - 2k) x + 5 pX {p + k) v , (4) 
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completes the description of the polarization tensor. These formulas hold also in a 
background field, provided the corresponding expressions for the propagators are 
used. We take a homogeneous magnetic background field in the representation 
where it is parallel to the third axes is both, configuration and color spaces, 



B S^ 3 . 



(5) 



After turning the momentum and the gluon field in color space into the charged 
basis, 



Pn = B^ a p a with B 



/iO: 





(6) 



we have a color neutral and a color charged field, both of spin 1. The charged field 
will occupy in the presence of the background field Landau levels and we expand 
with respect to the corresponding eigenfunctions. In this space, the momentum 
p^ of the charged gluon is an operator, whose components fulfill the commutation 
relation 



\Pa,Pp] = iF a/3 = i 



\ 









\ 

1 
1 
1/ 



B. 



(7) 



a/3 



In this expression we write B instead gB. In what follows, where it will be not 
misleading we also put B — 1, for short. 

A basis in this space is given by the vectors | n, a)^, eq. (30) in [IB] . The tree 
level energies of these states are 

l\ + ll + B(2n + 1 + 2a) 



in 



0,1, 



a 



±1), 



(8) 



where I3 and 1^ are the momenta in parallel to the background field and imaginary 
time, respectively. The lowest state is n = with o = — 1 and it is called 
tachyonic for allowing for negative energy, E n < 0. In the charged basis, the 
lowest {tachyonic) state is 



t) a =\ 0,-l) Q 



/ 1 \ 




V / 



0), 



(9) 



where | 0) is the lowest Landau level which is annihilated by the operator a in 

( ia) \ 

Pa = 



V 



-ia 
h 
h 



(10) 
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and we note 

Pf,\t)»=Pl\t)a = 0. (11) 

Below, we will also use the notations I 2 = I 2 + l\ and h 2 = p\ + p\ = aa) + a' a. 

Now, we continue with using the representation of the polarization tensor as 
given by eq. (51) in [13]. It results from the proper time representation of the 
propagators, 

oo oo 

G{p -k)= fds e~ s{p - k) \ G{k) = Jdt e~ tk \ (12) 



and integration over k in eq. (j2J). First we consider the case T = and drop 
the sum over iV in the equation mentioned. The representation is in terms of a 
parametric integral over the parameters s and t, 

n AV = f ds f dt e( s , t) (y: m% + Mf A h , j + n^ po1 (13) 

with 

e(M)= exp( ^V . (14) 

Here the following notations are used: 

st 

H = I 2 + m(s,t)h 2 , 

s + 1 

m[s,t) = s + - In — , 
2 ^ + 

H ± = t + sinh(s)e ±s , (15) 

which are equivalent to eqs. (23-26) in [13]. The sum over i, j in ( TT3l) follows the 
subdivision introduced in [13] and the functions M^{, are given by eq. (53) in 

m. 

Now we take the tachyonic projection of IIaa', eq. ( 1X3]) . In doing so we note 
especially h 2 = 1 (for B = 1) and the function simplifies, 

exof — ——I 2 — 

For the projection of the functions M A '|, we use representation (55) in [13]. We 
collect them in the appendix, together with the contributions from the tadpole 
graphs. At this place we mention that under the tachyonic projection we get 
directly a representation suitable for further calculations. This is in opposite to 
the general case in [13j, where we were forced to integrate by parts in a special 
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way in order to account for the transversal structure of the polarization tensor 
already on the level of the parametric integrals. 

From eq. ( 1571) in the appendix we get the non-regularized expression. Adding 
a dimensional regularization (e > with the renormalization mass scale /i) and 
restoring the dependence on the magnetic field B we get 

(t | n | t) = 




Using variables 

s = qu, t — q(l — u), (18) 

it is easily seen that the linear ultraviolet divergencies, which manifests them- 
selves in the behavior for q — > 0, cancel (this is a consequence of the gauge 
transformation properties of the polarization tensor [17]). The logarithmic ones 
remain and these take the form 




<*in|t) = 77Z^ - ^ -(l 2 -B)+0(l) . (19) 



(4tt. 



This is just the same as in case with zero magnetic background field (keeping, 
however, the basis of Landau level states). Hence we can use the standard renor- 
malization prescription as to subtract the polarization tensor with zero field. 
Specifically we define 



(t | n | t) div = 




0(1). (20) 



The factor 1/(1 + q) was introduced for convenience. This is an alternative to an 
auxiliary gluon mass. The arbitrariness brought in by this factor corresponds to 
a finite renormalization or a change in the parameter [i. 
Now we define the renormalized polarization tensor by 

(t | n | t) ren = (t | n 1 1) - (t | n | t) div k=0 . (21) 

Using ( TT7|) . (|2U|) and the variables (TTH|) . it takes the form 

B 



(t 1 n 1 ty 



POO 


dq { 


Jo 


V I 


_Q_ 


(— 


//_ 





(4tt) 2 
1 

du 



2 2 
q sinh(g) 
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This expression is ultraviolet finite. However, it has still a divergence resulting 
from the tachyonic mode in the loop integration. This divergence manifests itself 
in the s-integration for s — > oo. To understand this divergence one needs to 
remember that we already performed the Wick rotation. The initial expression 
needs to be written in Minkowski space with the corresponding '+ie '-prescriptions 
in the propagators. Starting from there, one needs to identify a contribution 
involving the poles (in momentum plane) from the tachyonic mode in the loop 
integration (tachyonic part). That contribution may be anti-Wick rotated (in 
order to keep convergence) which is q — > —iq in the parametric representation. 
The remaining part may be Wick rotated as usual (q — > iq). We mention that 
the identification of the tachyonic part is not unique which, however, does not 
influence the final result. In representation ( 122]) . once we already Wick rotated 
the whole expression, we can identify the tachyonic part and anti-Wick rotate it 
twice, substituting q — > qe ln . 

The tachyonic part is just that which, after Wick rotation is exponentially 
divergent for s — > oo. In order to identify this part we represent fl22|) in the form 

(t | II | t) ren = A 1 + A 2 + B 1 + B 2 , (23) 

where we defined from the double integral 

1 ralsdt te s 2 

Al -J^fl WtTVt 1 ~t + Y ( } 

and, from the single integral, 

Al * r d sm\. 2)e ,. (25) 



(4tt) 2 Jo q \ B 

These are the tachyonic parts, and B\ and B 2 collect the remaining contributions, 
see below. 

We mention that Ai has no ultraviolet divergence and we could put e = in 
( |24|) . Now we perform in these two expressions the anti-Wick rotation. We get 

_ st I 2 

1 rdsdt te~ s 2 e s + tB 

1 = J^yJ 7+t J+t t-i-io' ( } 

where the '— iO '-prescription results from turning the integration path. It makes 
Ai having an imaginary part. Now Ai, eq. ( |26|) . is a finite expression for I 2 < B 
and can be calculated directly, numerically for instance. 
The part A 2 , after anti-Wick rotation, becomes 

A2=7 ^r d s(fi3f\\. 2)e -, (27) 



(4tt) 2 7 q \ B 
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and can be calculated explicitly, 



A 2 = AT + Af\ 



with 



and 



.4: 



div 



B 



(47T) 



- + ln t " 



B 



(28) 
(29) 



AT 



-2ni, (30) 
where we made the subdivision according to the ultraviolet divergence in A 2 . 
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Figure 1: The function f im in eq. ( l38l) as function of its argument. 
Finally we have to calculate B\ and B 2 . We represent them in the form 



B l = BT + B{ 



div 



with 



B\ c 



B r dsdt j 
(4vr) 2 J 7+t { 



4 A s + te 

— + 4 

u_ s + t 



2s te~ s I 2 
4- 



s + tB 



g s+t B ' 



and 



div 



B r dsdt ( [i 2 q 
(An) 2 J 7+t \ W / 



s + t 



+ 4 



s(s-t) + (s 2 + t 2 + st)£ 
(s + t) 3 (l + s + t) 

s(s -t) + (s 2 + t 2 + st) 1 ^ 



{s + t) 3 (l + s + t) 



(31) 



(32) 



(33) 
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Figure 2: The function / re in eq. (138j) as function of its argument. 

for the part resulting from the double integral and 

B 2 = B l 2 cn + Bf\ (34) 



with 



(An) 2 Jo q \ q sinh(g) q 1 + q J 

2 - 54 (Jf (35) 



and 



?div B r d i (» <i\ ( - 4 



for the contribution from the integral over q. The divergent contributions collect 
into (EQJ), 

(t | n | t) div = Af v + £f v + Kf v , (37) 

and undergo the renormalization. 

In this way, we get for the renormalized polarization tensor 

Re((t\U\t)n = JL/ re ^ =Re(A 1 ) + Br + Br, 

Im((t | n I t)™) ee f^(^=Im(A 1 ) + Im(Ar). (38) 

The functions / re and /j m are shown in Fig.l. 
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3 Magnetic mass at high temperature 



The above formulas were obtained for zero temperature, T = 0. These allow 
for immediately writing down the corresponding results for high temperature, 
T — > oo. To do that, first one has to return to eq. (j5J). At finite temperature T, 
in that formula, the integration over the momentum fc 4 must be substituted by 
the Matsubara sum, 

°°fH T £ (39) 

-oo l=~co 

in the otherwise unchanged expression. From here, the limit of high tempera- 
ture, T — y oo, is obtained by taking the contribution of the zeroth Matsubara 
frequency, / = 0. This is, of course, the known dimensional reduction to a the- 
ory without temperature in three dimensions. In the representation in terms of 
parametric integrals, eq. (TTT|) . the integrand is simply 'missing' a factor \/4irq. 
Hence, the tachyonic projection of the polarization tensor (|T7|) takes the form 

(t | n | t) 

T->oo 




sinh(g) 




4cosh(g) }, (40) 



where now I 2 = I 2 . In this expression, we removed the regularization. This 
is possible since the ( |4"0]) does not contain ultraviolet divergencies. The former 
linear ones cancel as before (in four dimensions) and the former logarithmic ones 
do not appear due to the dimensionality. 

In the following we focus on the magnetic mass of the charged gluon, 

m 2 ch =< t\U(B,T,p 4 = 0,p 3 0)|t > . (41) 

It is given by eq. ( 14"0]) with I 2 = in the tachyonic projection. So it remains 
to treat the tachyonic mode in ( 140]) . Because of I 2 = this is here simpler 
than in the preceding section. The tachyonic contribution results from A2 (125]) 
only. Accounting for the changes dimensionality it can be calculated easily and 
it delivers the imaginary part, 



T / B 

Im((t I II I t)^) = -2iyfr—^, (42) 

The real part comes completely from the i? 2 -part. It is also much simpler than 
in the preceding section and it reads 



Re((t I II I t) 



T->ooJ 



1 4qe~ qu 2 2 _ 

(47r) 3 / 2 Jo y/q Jo fi 2 . q sinh(g) 

- 3 - 26 (4^' (43) 
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The last line is a result of numerical integration. 

The main result of this calculation is that the projection of the polarization 
tensor is proportional to Ta/5 allowing for the conclusion that there is a field 
depentent magnetic mass. Further, it is important that the real part has negative 
sign. Hence it follows that the radiation corrections act to stabilize the spectrum 
of charged gluons at high temperature. 

4 Conclusions 

In the above sections we calculated the tachyonic projection of the polarization 
tensor of color charged gluon in a chromomagnetic background field. We used the 
general formulas derived earlier in [13] . The expression for the polarization tensor 
is given on the level of the Feynman rules by eqs. (jSJ) and (JHJ). We account for the 
magnetic background field by the known representation in terms of parametric 
integrals after Schwinger. From simply inserting these formulas one comes to 
an expression where, already in QED, the transversal structure cannot be seen 
directly. Instead, one needs to integrate by parts in the parametric integrals to 
make this structure manifest in the integrands. A check for this procedure is 
the cancelation of the linear ultraviolet divergences. In the non-Abelian case the 
picture is more involved since the polarization tensor is not transversal as shown 
in [T7]. Nevertheless, in [13] a way was found how to perform the necessary 
partial integrations in order to make the corresponding structures, formulated in 
terms of form factors, manifest. 

In this way, in JT3] , the general structure of the polarization tensor in terms of 
the parametric integrals was found. However, as it turned out, this structure does 
not hold in general. For the projection on the tachyonic state, which we consider 
in the present paper, the integration by parts does not work. Doing it formally, 
artificial divergences are introduced which manifest themselves as factors /i + in 
the denominators which in fact should never appear. 

For this reason we made in the present paper a direct projection of the polar- 
ization tensor instead integration by parts. As a result we obtained the represen- 
tation (fTTj) . Here the linear ultraviolet divergence cancels between the different 
contributing graphs (the basis graph, the ghosts and the tadpole graphs). We 
consider this a check for having the correct properties with respect to the gauge 
symmetry. 

The remaining calculation involves the anti-Wick rotation and the renormal- 
ization according to the known rules. The result is quite simple and numerical 
examples are shown in Figs. 1 and 2. 

On the base of these calculations, we found the charged gluon magnetic mass 
in the high temperature limit. The derived result in eq. (l43p is signalling that 
radiation corrections act to stabilize the tree-level spectrum (jSJ). Really, if one 
considers the pole of the Schwinger- Dyson operator equation, taken in the ground 
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state, < t\D\t > _1 =< t\l 2 - gB-U(T, B,p 4 = 0,p 3 ->■ 0)|t >, then the positivity 
of the "effective mass squared" rr^^ = —gB + 3.26 ^f^ r at high temperature 
follows. This result fills in the gap existed in calculations of the screening parame- 
ters for charged gluons. Due to spectrum stability, the stabilization of the Abelian 
chromomagnetic fields of the type in eq.© at high temperature is guaranteed. 

These results can find a number of applications in different fields. At finite 
temperature, these are the early Universe and heavy ion collisions where differ- 
ent types of magnetic fields are expected to be generated. At zero temperature, 
the obtained results may serve as starting point for the calculation of various 
characteristics of charged vector particles with gyromagnetic ratio 7 = 2. An ex- 
ample is the p-meson electrodynamics in strong magnetic fields, where interesting 
properties of the vacuum were discussed in the literature (see, for example, [18J, 
[19] . (20] and references therein). Another obvious candidate for application is a 
M^-boson in magnetic fields. In the two latter cases, one needs to account for the 
mass term of the particle, which can be incorporated in the developed formalism 
easily. 
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Appendix 

In this appendix we perform the projection of the functions M l x \,, eq. ( fT3l) . 
using the expressions given in eq. (55) in [13]. We repeat that formulas here 
putting u = and stripping off a factor Z since we do not need it in the current 
calculation. In [13J it was introduced for use in integrating by parts which we 
do not need to do here. Furthermore, we drop all contributions containing p\ in 
view of eq. ffTTj) and obtain 



M- 



rll 

AA' 




M- 



-12 

AA' 



M 



■21 _ 

AA' _ 








M- 



r22 

AA' 



4S XX , (pT T p) + 4T A T V (pSp) + AST T 2iF , 
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r2:-! , r? n 2 — ' 2 «F \ / 2«F 



M$ = M&f = -2tS w (PprP) ~ IMF D 



M% = Mfg = -2E M (p(±) p) + ( A ^) xy ■ ( 44 ) 

Here notations from [13] are used, which we rewrite somehow, 

2iF 5u f-1 1 \ / 1 1 \ _ 

— f- = — ^ + — + — \iF+\ — + — )5 ± , 

L) 1 S + t fj, + J \fi_ fi + J 

A s . f 1 — a 1 — a\ ( 1 — a 1 — a A . 

— = 5ii + \i>F+\ 1 ]6±, 

D s + t " V 4//_ 4a/i +y / V 4/W- 4a/i+ / 

y /i_ a/x + y y //_ y 

r = dn + pit-fcW+f&.fcV, ( 45) 

and we introduced the new notations, 

a = e" 2s , /i± = i + sinh(s) e ±s , (46) 
and we remind eq. (20) in [13], 

E = S\\- sinh(2s)iF + cosh(2s)<5 ± . (47) 

In order to calculate the contributions quadratic in the momenta we note the 
formula 

(p{fiS\\ + aiF + u5±)pj = [il 2 - o + v, (48) 

which holds for any fi, a, v in the tachyonic state (h 2 = 1). From here we 
calculate 



( 2tF \ 


I 2 1 




S + t /i_ 


® T >) 


S 72 1 

= r + - 

s + t : 


(pT T p) 


= / 2 + ^i, 




(pSp) 


^2 + 



2/x. 



(49) 

Further we need the tachyonic projection, 

(t | fi8\\ + aiF + v5 L | t) = -a + v , (50) 
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which follows directly from eqs. (29) and (32) in [13]. Now, using (149]) and ([50 

1 + cosh(2s) 



we can calculate the averages in (Tj 


HD, 




<* 


\M n 


*> = 


(t 1 M 131 1 t) + (t 1 

\ 1 1 / ' \ 1 


M Sl1 


t) = 

1 


(t | M 132 | t) + (t | 


M 3al 


t) = 




| M 22 


t) = 


(t | M 231 | t) + (t | 


M 3x2 


t) = 


(t | M 232 | t) + (t | 


M 3 2 2 


t) = 



a\i- 

-— ' 2 , 

s + t 
4s Z 2 

s + t a v ' 

In these formulas some simplifications occurred. Now we add these contributions, 



2a 
li- 



lt I 1 1) = A_ 4 ^±l /2 

1 ^ 17 s + t 



(52) 



where, again, a number of contributions canceled. It remains to consider the 
contribution from M 33 and from the ghosts (last line in eq. (55) in [13J). These 
combine into derivatives, see eq. (90) in [T3], which allows for carrying out one 
of the parameter integrations. Using 

6(s = 0, t) = — , 6(s, t = 0) = , - - - , (53) 



t 2 ' 



s sinh(s) 



we get in the projection 



/ dsdt (t | (M 33 + M gh ) 9(s,t) \t) = [—(- + — V— ) . 
J v ' J q \q asmh(q) J 



(54) 



Finally, we need the contributions from the tadpoles, ([3]), which take the form 

(55) 

Together, ( 154"]) and f l55|) collect into 



(t I n tp It) = I — (—- 2 + cosh ( 2 ^ g ) ~ 3 sinh ( 2 g) \ 
q \ q sinh(g) J 



(t | n 33 + n gh + n tp 1 1) 



1 rdqf-2 1 + sinh(2s) 



(4tt) 



Collecting from fl52|) and fl56|) we come to 



(t\U\t) 



1 



(4tt) 2 
+ 



c?s dt 
s + t 
-2 



^ — -2- 
q \ q 



4 s + 1 
^4 SL 

//_ sinh(s) 



sinh(s) 



(56) 



__s^,2_ 

e s+t 



q \ q sinh(s^ 
14 



— 4 cosh(g) 



(57) 



which represents the non-renormalized polarization tensor of the charged gluon 
in the projection to the lowest, the tachyonic state. 
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